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Abstract. We have investigated the formation of helium droplets in two physical situations. In the first one, 
droplets are atomised from superfiuid or normal liquid by a fast helium vapour flow. In the second, droplets 
of normal liquid are formed inside porous glasses during the process of helium condensation. The context, 
aims, and results of these experiments are reviewed, with focus on the specificity of light scattering by 
helium. In particular, we discuss how, for different reasons, the closeness to unity of the index of refraction 
of helium allows in both cases to minimise the problem of multiple scattering and obtain results which it 
would not be possible to get using other fiuids. 

PACS. 67.25. bf Normal phase of 4He : Transport, hydrodynamics - 67.25.bh Films and restricted ge- 
ometries - 47.80.Jk Flow visualization and imaging - 47.55.-t Multiphase and stratified flows - 67.25.dg 
Superfluid phase of 4He : Transport, hydrodynamics, and superfiow 

1 Introduction eating objects inside turbid media has been shown to be 

possible by measuring the statical or dynamic scattered 
Light scattering is a widely used tool to investigate prob- intensity. In particular, the so-called dynamic wave spectroscopy [1 
lems in soft matter physics. In many instances, multiple [3] enables the visualisation of flows[l] or the location of ab- 
scattering turns out to be a problem. During the last sorbing objects |5J inside multiple-scattering liquids. How- 
twenty years, significant progresses in dealing with and ever, imaging under multiple scattering conditions is chal- 
even exploiting multiple scattering have been made. Lo- longing and information provided by such techniques re- 
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tering. This explains why a great deal of effort has been This implies that the scattering cross-section for scatterers 

devoted to the suppression or the minimisation of multiple much smaller than the light wavelength is very small, so 

scattering in dense soft-matter systems. A first example that multiple scattering barely occurs. The situation for 

is the study of the glass transition in colloidal suspen- larger scatterers is more subtle. In the regime of geomet- 

sions of silica and copolymer spherical particles[6j. Here, rical optics, the scattering cross section is approximately 

the solvent index of refraction is matched to that of the twice the geometrical cross-section (diffraction and reflec- 

copolymer particles, so that scattering only arises from tion/refraction contributing nearly equally[Tl|), so that, 

the silica particles. Their volume fraction is small enough at a given volume fraction, the scattering mean free path 

(2%) to study their dynamics using quasi-elastic light sin- decreases with the particle diameter (figure[l]), and reaches 

gle scattering (the small amount of multiple scattering a minimum for a particle radius of about 10 /im. Even for 

being eliminated by two-colour cross-correlation [7j). A moderate volume fractions (10"'^), it can become smaller 

second example is the study of phase separation of binary than a light path length of 1 cm, giving rise to multiple 

mixtures in the presence of a disordered environment, pro- scattering. However, the small refractive index of helium 

vided by confining the mixture inside a porous medium[51 implies that the scattering is nearly forward in this case. 

[51ll0|. Near the critical point, the density fluctuations give Figure [1] shows that, for droplets of liquid helium, the 

rise to multiple scattering which is minimised by careful transport mean free path /* is then more than two orders 

choice of the refractive index of the fluids with respect of magnitude larger than the scattering mean free path 

to that of the substrate (e.g. isobutyric acid and water in Igca- As we will show, this peculiar situation makes imag- 

silica aerogels |10j). ing possible even through samples much larger than Igca- 

Such an index matching does not occur when one con- In this paper, we will describe two physical situations 

siders a diphasic pure fluid. In this case, the optical con- where we take full advantage of the low refractive index of 

trast between the vapour and the liquid phases is solely helium. The first one is the study of atomisation of nor- 

fixed by the ratio of their densities and cannot be tuned, mal or superfluid liquid helium by a fast vapour flow. The 

This makes the optical study of diphasic systems prone to droplets created through this process range from microns 

the occurrence of multiple scattering. From this point of to tens of microns, so that geometrical optics applies. Even 

view, helium is a remarkable exception. Due to its small when the droplets density becomes large, observation at 

size and closed shell electronic structure, its polarisability appropriate angles allows to determine their spatial distri- 

is very low and its index of refraction close to 1. For liquid bution, and to analyse the influence of the physical charac- 

helium below 2 K (density of 145 kg/m'^), the dielectric teristics of helium on the atomisation process. The second 

permitivity is 1.057[TT], and the index of refraction 1.028. situation is the condensation and evaporation of helium 
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Fig. 1. Scattering and transport mean free paths of light Isca 
and r (in mm), for uncorrelated helium and water spherical 
droplets of radius a (index contrast 1.025 and 1.33), at a vol- 
ume fraction of 5.10~^. The wavelength is 632 nm. 

inside porous silica glasses of different topologies, Vycor 
and silica aerogels. Both processes involve fluctuations of 
the helium density on a subniicronic scale, which can be 
probed by spatially resolved light single scattering. This 
provides information on the size of the microscopic het- 
erogeneities and their macroscopic distribution, which can 
both be compared to the predictions of the existing models 
for the condensation and evaporation processes. 

2 The case of large scatterers : Atomisation 
of liquid helium 

2.1 The Cryoloop experiment 

Blowing a fast enough gas stream parallel to the free sur- 
face of a liquid results in the liquid atomisation into a 
spray of droplets. The Cryoloop experiment aims at study- 
ing this process for diphasic helium in a horizontal pipe. 



Table 1. Comparison of the liquid and vapour densities {pL, 
pv) and viscosities {tjl, rjv), and of the liquid- vapour interfa- 
cial energy (a) for the diphasic helium and water-air systems; 
data for helium from refs. [T3| and (for a) !14!. 



fluids 



PL 



pv 



VL 



T]v 



(kg/m^) (kg/m^) (^Pa.s) (^Pa.s) (mj/m^) 



He 1.8 K 


145 


0.48 


1.3 


0.43 


0.316 


He 2.0 K 


145 


0.79 


1.5 


0.52 


0.301 


He 2.2 K 


146 


1.25 


2.6 


0.57 


0.284 


He 2.6 K 


144 


2.55 


3.3 


0.69 


0.253 


Water/air 


1000 


1.2 


1000 


17 


70 



The liquid can be either normal or superfluid, depending 
on the temperature, and the gas is helium vapour at the 
saturated vapour pressure. The interest of this experiment 
is two-fold. First, diphasic flows of helium can be used for 
refrigeration purposes, an example (below the atomisation 
threshold) being the cooling of the superconducting mag- 
nets of the Large Hadron Collider, and it is important 
to assay the conditions and consequences of atomisation. 
Second, the physical parameters relevant for atomisation 
are very different for diphasic helium and for the widely 
studied air- water system (see table [ij , which offers an op- 
portunity to experimentally test the influence of these pa- 
rameters on the atomisation process. 

The Cryoloop facility is installed at the Service des 
Basses Temperatures at CEA-Grenoble. A specially de- 
veloped refrigerator [TSj with large cooling power (400 W 
at 1.8 K) is used to deliver a flow of liquid of up to approx- 
imately 20 g/s (i.e. more than 400 litres of liquid per hour) 
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to one end of a nearly horizontal, 10 ni long and 40 mm 
inner diameter, pipe. This liquid is partly evaporated by a 
heater in order to obtain downstream an essentially strat- 
ified diphasic flow in which the bulk liquid moves along 
the bottom of the pipe, and the vapour moves above the 
liquid. At the end of the pipe, the remaining liquid is evap- 
orated, and the resulting vapour is pumped through cold 
compressors followed by a room temperature pump, which 
sets the pressure, hence the temperature, inside the pipe. 
At a given applied power (i.e. a given mass flow rate of 
vapour), the vapour velocity depends on the cross-section 

of the pipe accessible to the vapour and on its density, tt « tt i- • ^ ^ 1 1 ^ ■ x- • a i- 

'^ "^ ^ -^ 1? ig. 2. Helium mist generated by atomisation in a 4 cm diame- 

In our experiments, the former is essentially the whole ^^^ p-p^ (iog_gray scale). The inner walls of the pipe correspond 

cross-section of the pipe, so that the ultimate maximal to the white bars. The mist is visualised using illumination by a 

(mean) vapour velocity decreases with increasing tem- laser sheet at 15° from the CCD direction of observation (left), 

perature (since the saturated vapour pressure, hence the and by a vertical laser (right). The liquid- vapour interface is 

vapour density, increases with temperature). In order to visible as a brighter region in the bottom of the pipe, 
obtain atomisation, a part of the total flow has to be kept 

liquid, which makes the practical maximal velocity to de- 
ports through the cryostat containing the pipe allow to il- 

pend on the amount of liquid left. Typically, for a liquid 

luminate and to detect the spray using room-temperature 

height of 3.5 mm (corresponding to a wetted fraction of 

optical components 16J. The spray distribution across the 

the pipe of 20%), the maximal velocity decreases from 

whole cross-section of the pipe can be observed by illumi- 

18 m/s at 1.8 K to 6 m/s at 2.5 K. 

nating the pipe with a laser sheet propagating horizontally 

perpendicular to the pipe. A CCD located opposite to the 
2.2 Visualisation of atomisation 

laser source with respect to the pipe is used to image the 

In order to observe the droplets created by the atomisation illuminated section under a 15° angle, as shown in figure[21 

process, the stainless steel pipe is replaced along a length With the same CCD, we also observe the light scattered 

of 13 cm by a Pyrex tube of identical inner diameter. This at 90° from a laser beam propagating vertically along the 

optical section is located close to the end of the pipe to pipe diameter, which allows to compare scattering at two 

ensure that the observed spray is fully developed. View- angles along this particular light path. 
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Fig. 3. Angular dependence of the intensity scattered by the 
mist, as measured on a screen perpendicular to the incident 
laser beam. The unscattered beam passes through the hole at 
the center of the picture. Measuring its intensity gives the scat- 
tering mean free path, averaged over the horizontal diameter of 
the pipe. White spots correspond to reflections on the windows. 
The scattered intensity does not depend on the azimuthal an- 
gle <j!>, which means that the droplets are not elongated by the 
gas flow. Averaging this intensity over (j) gives the thick curve. 
Above 40 mrad, it follows the prediction of geometrical optics, 
showing that the droplet size is larger than one micron. The 
rise at smaller angles is due to diffraction and compares well to 
the Mie prediction for an exponential distribution of diameters 
with characteristic diameter 20 ^m. 

A second optical access opens to a stainless steel por- 
tion of the pipe through small fiat windows, centred on the 
pipe horizontal median plane. Shining a horizontal laser 
beam through these windows and visualising the scattered 
intensity on a screen located on the other side of the pipe 
allows to determine the angular dependence of the scat- 
tered intensity in a cone of about 10° around the forward 



direction. This enables to test the scattering regime, i.e. 
the size of the droplets, on the median plane of the pipe. 
Figure [3] shows that the scattering is strongly peaked in 
the forward direction. The measured angular dependence 
is that predicted by geometrical optics for scattering an- 
gles larger than 40 mrad, with a steep rise at smaller an- 
gles, associated with diffraction. The observation of refrac- 
tion shows that the average diameter of droplets is larger 
than one micron. This is confirmed by both the angu- 
lar range over which diffraction occurs and measurements 
with a Phase Doppler Particle Analyser (PDPA), which 
show that the average diameter is a few tens of microns, 
for all hydrodynamic conditions. 

Being in the regime of geometrical optics implies that 
the scattered intensity (outside the diffraction peak) is 
directly proportional to the interfacial area S, per unit 
total volume, of the droplets. This also holds for the scat- 
tering mean free path, since the total scattering cross sec- 
tion is twice the geometrical one (more precisely, there are 
oscillations with size, due to interference effects, of both 
quantities, but they are washed out by the average over 
the droplets size distribution). By comparing, for given 
conditions, the image of the mist with the attenuation of 
the horizontal beam used for measuring the angular de- 
pendence, we can determine the coefficient between the 
brightness of the images and S, and thus, extract from 
any image the spatial dependence of S (or, equivalently, 
of the local mean free path of light Isca—'^/^)- 

Figure [4] displays pictures separately recorded with 
laser sheet and vertical laser illumination for increasing 
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Fig. 4. Mist for velocities (averaged over the pipe section) of 
7, 10.5, 14, 17.2 and 17.6 ni/s, at a temperature of 1.8 K and 
a liquid level corresponding to the bottom 20% of the pipe 
perimeter. The images are averages of several tens of pictures 
with 30 ms exposure time (10 ms for the two pictures above 
17 m/s). A log-gray scale is used both for sheet and vertical 
laser illumination. 

vapour velocities, at constant bulk liquid level and vapour 
temperature. The obvious strong increase of atomisation 
with vapour velocity is quantified in figure O which shows 
the dependence of the interfacial area on the height z in- 
side the pipe, as deduced from both modes of illumination. 
The curves show that the mist is stratified over a char- 
acteristic height Hchar which increases with the vapour 
velocity. This stratification results from the competition 
between the diffusion of the droplets due to the vapour 
turbulence and their gravitational fall through the viscous 
vapour [T7]. The decrease of stratification with increasing 
vapour velocity can then be ascribed to the increase of the 




40 z (mm) 



Fig. 5. Mist stratification ; Profiles of interfacial density de- 
duced from figure|3]for sheet (15°, top) and vertical laser (90°, 
bottom) illumination. The left scales gives the corresponding 
ratio of the pipe diameter D to the local mean free path Isca- 
The right scale is l/hca = S/2. The increase of signal beyond 
40 nun is due to scattering from the pipe walls. The peaks 
above 20 mm for laser illumination arise from similar scatter- 
ing of the beam reflected by the bottom of the pipe. 

droplets turbulent diffusivity or to a decrease of their set- 
tling velocity due to a smaller size. Beyond this common 
behaviour, a striking difference between the two modes of 
illumination is that, in the bottom part of the tube (z < 15 
mm), the interfacial intensity deduced from sheet illumi- 
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nation saturates while that deduced from the vertical laser to multiplications if !?'(«) and ^^^'{u) are expanded on 
still grows. This discrepancy occurs for values of the mean Legendre polynomials : 
free path smaller than the pipe diameter (which is close 



to the distance travelled by the laser sheet, except very 
near the bottom of the pipe). This suggests that multi- 
ple scattering is involved. The question then arises of how 
the path of the vertical laser can be imaged in a region 



'P{u)^Y.^2l + l)viPi{u) (1) 

I 

OQ 

W^P\u)^Y.^2l + l){viYPi{u) (2) 



1th 



where the mean free path is as small as a tenth of the pipe ^' ^/^ / ^{u)fi{u)du, [^6) 

diameter. A quahtative answer is that, because the trans- ^j^j^ ^^^^g ^^ compute ^'^p\u) for any given initial 

port mean free path for helium droplets is so much larger distribution ^{J^ A particular case is that of the phase 

than the scattering mean free path, the hght scattered at f^^ctions M^^iu, v) corresponding to coefficients vi given 

a large angle involves only one scattering event at (ap- i^ . 
proximately) this angle, between two series of scattering 



events at small angle. However, it is not obvious whether 
this mechanism is more likely than a series of small-angle 
scattering events. In fact, as we now discuss, this depends 
not only on the ratio of the transport mean free path I* 
to the scattering mean free path Isca, but also [TB] on the 
shape of the phase function, which describes the angular 
dependence of the scattered intensity for one scattering 
event. 



yi = {vf (4) 

where 1 < a < 2, and v =< u >, the average of cos 6* over 
the distribution 'F{cos9). Note that 1 — v — Isca/l*, and 
that strong forward scattering corresponds to v close to 1. 
These phase functions have a functional form stable under 
the convolution process. Special cases are a=l and a=2, 
corresponding respectively to the Henyey-Greenstein (H- 
G) phase function [5^ : 



2.3 The multiple scattering of light by helium droplets ^ (^, ,,\ ' "" > (k\ 

The scattered intensity by a spherical droplet only de- ^nd to a nearly Gaussian distribution (the true Gaussian 
pends on the angle 9 between the incident and the scat- distribution corresponds to i^e = w^(^+i)/2). In both cases. 



tered directions, and is described by the phase function i Qne should stress that this represents the angular distri- 

!f (u =cos 61), the probability distribution of 9 after one bution of the light, when summed over all possible positions 

scattering event, multiplied by An. The probability distri- after the p scattering events, i.e. the exact equation [2] does not 

bution 'P^P' (u) after p scattering events is obtained by sue- give any information on the spatial distribution of intensity, in 

cessive convolutions of <F{u). These convolutions amount contrast to the full radiative transfer equation[19]. 
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the angular dependence after p scattering events is given 
by the same functional form, with Vp — v^. 

For small angular width and scattered angles, the H-G 
function is closely similar to a 2D Lorentzian distribution 
around the forward direction. For a given half width at 
half maximum (HWHM) Ohw, the large "wings" of such 
a distribution result in an r.m.s. value much larger than 
9hw, Orms =< 0^ >^/'= {2lsca/l*y/^ of Order (^w)'/' 
This is in contrast with the Gaussian phase function of 
same 9rms for which the HWHM is of order Orms- As first 
pointed out by Van de Hulst [21] . this implies that the H-G 
phase function behaves very differently from the Gaussian 
phase function of same 9rms under multiple scattering. In 
both cases, Vp = v^ implies that 9rms scales with (p)^'^, 
but the characteristic width for the H-G case increases 
linearly with p, i.e. much faster than in the Gaussian case 
(p^/^). As shown by figure[S]Ja), for a ~ 3/2, intermediate 
between the H-G and the Gaussian cases, 9hw scales as 
p^l'^ ^ suggesting that the general behaviour is of the form 
9hw k P^' ", although we have not proven it analytically. 

Equation [5] also gives the behaviour at large angles. 
Figures [nijb) , [DJc), [Sfd) show how the intensity multiply 
scattered at various angles 9 varies with the average order 
of scattering < p >, for a value of /*// = 1/(1 — v) = 
200, a ratio close to that for helium, and respective val- 
ues of a of 1, 1.5, 2. Here, < p > is the average num- 
ber of scattering events, and the corresponding multiply 
scattered phase function is obtained by summing equa- 
tion [5] over all p, with a weight given by the Poisson 
distribution associated with < p >[TF. In all cases, the 
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Fig. 6. (a) Half width at half maximum of the angular distri- 
bution of intensity, after p scattering events, for three single- 
scattering phase functions of same r.m.s. width < 9 > — 
2lsca/l*~0.01, differing by their values of a (see text). The 
straight lines correspond to a law p^^"'. The widening of the 
angular distribution increases when a increases, i.e. when the 
original phase function has extended "wings"; (b), (c), (d) : 
scattered intensity at different angles as a function of < p > 
(see text) for the different values of a. At large angles, the scat- 
tered intensity is initially linear in < p > for a=l and a— 1.5. 
This corresponds to a pseudo single-scattering regime. 

intensity is transferred from small to large angles until 
it is completely isotropized, but the way this isotropiza- 
tion takes place strongly depends on the angular shape 
of the original phase function, i.e. on a. In the Gaus- 
sian case, the intensity at large angles only comes from 
the accumulation of small single-scattering events : the 
larger the angle, the larger the number of events needed 
to give some signal at this angle. In contrast, for the 
H-G case, the signal at large angles (compared to the 
HWHM of the original phase function) initially grows lin- 
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Fig. 7. (a) Comparison of the Mie phase function for an expo- 
nential distribution of diameters with a characteristic diameter 
of 20 /im (at a wavelength 632 nm and for an index of refraction 
contrast of 1.025) to *'3/2(cos 0) for r/Zsca=250; (b) Intensity 
multiply scattered at various angles as a function of the aver- 
age order of scattering < p > for the Mie phase function. Note 
the existence of a linear regime up to 30 scattering events for 
a scattering angle of 90° . 

early with p, whatever the angle. This implies that it 
comes from one large angle single-scattering event sep- 
arating two sequences of small-angle scattering events, 
which nearly preserve the direction of propagation. Be- 
cause the large angle single-scattering event can occur at 
any scattering step, the intensity is indeed proportional to 
p, the total number of steps. The range of validity of this 
regime increases with the observation angle, and reaches 
about 200 scattering events at 90°. At this angle, a pseudo 
single-scattering behaviour can be observed over a dis- 
tance of order the transport mean path. The extension of 
the pseudo single-scattering regime depends on the shape 
of the phase function. Figure [6jd) shows that this regime 
still exists for a=1.5, but in a smaller range (< p > <10, 
for angles 6* >20°). 

What about the real case of helium ? Figure[7]Ja) shows 
the Mie phase function for a realistic distribution of he- 



lium droplets. The pronounced decay with angle is ap- 
proximately accounted for by classical refraction, below 
its maximal angle of deviation (w 20°), and by reflections 
above, the bump around 45° being due to the rainbow 
effect. The rise at small angles (inset) is due to diffrac- 
tion, or, more precisely, anomalous scattering[T2]. The Mie 
phase function is approximately represented by !f'3 /2(u) for 
I* /lsca—'250, so what we expect a behaviour similar to that 
of figure[6Ud). This can be checked by numerical computa- 
tion of the coefficients vi in this case and use of equation [5] 
to obtain the intensity multiply scattered at various angles 
as a function of the order of scattering < p >. Figure [7{b) 
confirms our expectation of a linear dependance in < p >, 
over a range which increases with the observation angle. 
In practice, we expect the pseudo single-scattering effect 
to hold up to 10 scattering events for an observation angle 
of 15° (case of the sheet illumination) and up to 30 scat- 
tering events for an observation angle of 90° (case of the 
vertical laser). 

In our experiment, we have tested these calculations up 
to < p >= 3 by comparing the intensity scattered from 
the pipe horizontal median plane to the average number 
of scattering events < p > deduced from the attenuation 
of the horizontal laser beam. Figure [5Ja) shows that the 
intensity scattered at « 4° as measured on the screen (fig- 
ure [3]) in not linear with < p >. Its dependence is well 
accounted for by our calculation. In contrast, figure [8jb) 
shows that the intensity scattered at 15° and 90° as mea- 
sured on the images using sheet and vertical laser illumina- 
tion, respectively, linearly depends on < p >, as expected. 
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Fig. 8. (a) Scattered intensity (arbitrary units) at 4° as a func- 
tion of the average number of scattering events < p > deduced 
from the attenuation of the beam. The incident laser beam is 
horizontal and probes the median plane of the tube. Different 
symbols correspond to different experimental conditions. The 
continuous curve is the predicted behaviour, (b) Scattered in- 
tensity at 15° and 90° on the pipe axis as a function of < p >. 
Data for the vertical laser beam are vertically shifted for clar- 
ity. The linear behaviour at < p > larger than 1 is specifically 
due to the strong forward scattering of helium. 

Based on figure [TJb), we would expect the sheet illu- 
mination to give the same interfacial area as the vertical 
laser for values of D/lsca of order 10. Figure [5] shows that 
this is not the case. This disagreement is related to the fact 
that our calculation predicts the total intensity scattered, 
but not where it does come from. Small scattering-angle 
events widen the incident beam, which spreads the signal 
on the CCD over a larger region than for simple scattering. 
This effect is clearly visible on the figure [5] in the bottom 
of the pipe, both with the vertical laser beam, and with 
an horizontal one. Thus, measuring the local interfacial 
area from the brightness of a pixel on the axis of the hori- 
zontal laser beam leads to underestimate it. As shown by 
figure ini^c), the larger the interfacial area (i.e. the closer 
to the bottom of the pipe), the larger the underestima- 
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Fig. 9. Comparison of the vertical profile of interfacial den- 
sity measured by different methods : laser sheet, and vertically 
scanned horizontal laser beam observed at 15°, vertical laser 
observed at 90°. (a) is a superposition of four images of the 
horizontal laser beam at four different heights. In the bottom 
part of the tube, where the droplets density is the largest, the 
beam is widened by multiple scattering. The same effect affects 
the vertical beam (b). Correct measurements of the interfacial 
area require to integrate the brightness over the beam width 
(box in (a)) rather than to take its maximal value on the axis 
on the beam. The results at 15° and 90° are then consistent 
(c). 

tion. The error made with the sheet illumination is less 
than with the horizontal laser beam, because the spread- 
ing of light outside a given pixel is partly compensated 
by the spreading from lower or larger heights. However, 
unlike with the sheet, a valid measurement of the interfa- 
cial area remains possible with the horizontal laser beam 
by integrating the optical signal over the laser whole ef- 
fective width, as shown in figure [9jc) (such an integration 
was also performed to obtain the points of figure [5{b)). In 
this way, the interfacial area deduced from the scattering 
at 15° is consistent with that deduced at 90°. 
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To summarise this section, we have shown that mean- 
ingful measurements of the intcrfacial intensity are possi- 
ble with helium, even in a regime of quite strong scatter- 
ing, provided that one observes at a large angle enough 
and one takes into account the widening of the beam due 
to multiple scattering. This advantage of helium was cru- 
cial to our optical measurements of developed atomisation. 

2.4 Atomisation threshold and stratification 

The combination of optics with the 400 W refrigerator ca- 
pabilities enabled us to study the dependence of the atom- 
isation process on three parameters : the level of liquid in 
the pipe, the vapour velocity, and the temperature, the 
latter setting the density of the vapour, and the state of 
the liquid (superfluid below 2.17 K, normal above). The 
results[22] will be described elsewhere[^. In particular, 
we have characterised two properties which determine the 
potentialities of the atomised flow for refrigeration, the ve- 
locity threshold for atomisation, and, above this threshold, 
the mist stratification. In the following, we will describe 
these properties and compare them to those for the water- 
air system, so as to pinpoint the specificities of helium. 

Let us first discuss the atomisation threshold, and its 
sensitivity to the superfluidity of the liquid. Figure [TU] 
shows how the the PDPA counting rate increases with 
the vapour velocity at constant temperature and total 
flow-rate, both below and above the superfluid transition 
(1.8 K and 2.52 K). For both temperatures, droplets are 
detected on the pipe axis above the minimal velocity used, 
about 4 m/s. The exponential decrease of the counting 
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Fig. 10. The vapour velocity dependence of the PDPA count 
rate gives an upper bound of 4 m/s for the atomisation thresh- 
old. The section probed by the PDPA (about 150 x 500 fim) 
is located on the pipe axis. The velocity is the average droplet 
velocity, as measured by the PDPA, which is about 20% larger 
than the average of the vapour velocity over the pipe section, 
computed from the vapour total flow rate. The total helium 
flow is 14.5 g/s at 1.8 K, and 16 g/s at 2.52 K, correspond- 
ing to a bulk liquid level at 4.5 m/s of approximately 35% the 
pipe perimeter at 1.8 K, and 24% at 2.52 K. Despite the larger 
vapour density at 2.52 K, the count rate is lower, due to the 
lower level, combined with the stronger stratification. 

rate at small velocities does not allow to extrapolate a 
velocity threshold from these measurements. The atomi- 
sation threshold was thus determined from pictures of the 
whole pipe section using sheet illumination. Droplets were 
detected above 3.4±0.4 m/s at 1.8 K, and 2.5±0.3 m/s 
at 2.52 K. These velocities are two to three times the 
Kelvin- Helmoltz threshold for wave formation [24j, Ukh = 
{iagpL/pyY^'^, where g is the gravity acceleration, and 
the other quantities are deflned in table[T] Ukh is 1.7 m/s 
at 1.8 K et 0.8 m/s at 2.52 K. They can also be compared 
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to the Ishii-Grolmes criterion[25J for atoniisation, which 
takes into account the effect of viscosity in a phcnomeno- 
logical way. It predicts a threshold of about 2.2 m/s at 
both temperatures, the smaller vapour density at 1.8 K 
being compensated for by the smaller viscosity. While it 
remarkably accounts for the threshold in the normal state, 
the Ishii-Grolmes criterion underestimates the threshold 
in the superfluid state, and, worse, predicts an atomisa- 
tion threshold below the threshold for wave formation, 
showing that it overestimates the influence of a vanishing 
viscosity. In any case, our experiment shows that the in- 
fluence of superfluidity on the atoniisation threshold is at 
most limited. The threshold is lower than for the water- 
air system (where it typically lies in the range 10-15 m/s, 
depending on the pipe diameter and the liquid flow rate 
[25l[26] ). but this can be mainly ascribed to the small den- 
sity and surface tension of liquid helium, the product of 
which is about 1/2000 of that for water. As this product 
enters Ukh with a power 1/4, we also understand why the 
difference between helium and water remains nevertheless 
limited. 

Above the atoniisation threshold, a stratifled mist is 
observed. Measurements at constant velocity show that 
the temperature (in the superfluid region) and liquid level 
do affect the degree of atomisation, but not the profile 
of stratification. This allows to study the effect of veloc- 
ity on stratification by comparing experiments at different 
temperatures and flow rates. 

Figure [TIT a') shows the stratification profile for three 
temperatures in the normal region, as deduced from the 
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Fig. 11. (a) Stratification for normal fluid, as measured using 
the vertical laser. The applied power is 360 W, and the total 
helium flow rate is 23 g/s. The liquid level covers from 22% to 
24% of the perimeter. The inverse of the slope of the straight 
lines is the stratification height Hchar- (b) Velocity dependence 
of H*f^^^—riv{1.8 K)Hchar/riv{T), which cancels the theoreti- 
cal dependence of Hchar on the vapour viscosity, for points at 
1.8 K (squares) and in the normal phase (circles). The vapour 
velocity Vv is measured by the PDPA on the pipe axis. The 
straight line corresponds to a V^ dependence. 

intensity scattered at 90° from the vertical laser, which 
we have shown in H2.3l to be the least sensitive to multiple 
scattering effects. The applied power being constant, the 
vapour velocity decreases with increasing vapour density, 
hence temperature. As in the superfluid phase, this implies 
a stronger stratiflcation. For both flgureslSJandfTTl the pro- 
file is approximately exponential, so that we can quantify 
the stratification by the characteristic length Hchar such 
that i:{z) oc exp{-z/ Hchar)- Figure [TT^b) shows Hchar, 
measured above 10 mm, for the profiles of figures [5] and 
lll[ as a function of the vapour velocity V^ measured on 
the pipe axis with the PDPA. Hchar is found to follow an 
approximate V^ dependence, similar to what we have ex- 
tracted from droplets concentration profiles measured for 
the water-air system[37]. 
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For droplets of diameter d small enough to follow the smaller than reported for the water-air interface. For ex- 
gas turbulence, Hchar ~ Dturb/Vfaii, the ratio of the tur- ample, for the latter system, Simmons and Hanratty |28| 
bulent diffusivity Dturb to the droplets fall velocity Vfaii oz report average diameters on the pipe axis decreasing from 
cP /rjy. For isotropic turbulence, Dturb is expected to scale 60 to 35 /im, for air velocities increasing from 30 to 50 
with the vapour velocity, so that iJc/iar ccrjyVy/cP. Hchar /vvt^/s (i-e- two to three times the atomisation threshold), 
vary approximately as Vy, suggesting that the droplets di- Here again, the difference probably stems from the very 

— 1/2 

ameter decreases with the vapour velocity as d ex K, , low surface tension of helium, 
independent on the superfluid nature of the liquid. Also, 

since the stratification, at a fixed vapour velocity, does not 3 -p^e case of Small scatterers : condensation 

depend on temperature, the droplet diameter should only 

and evaporation of helium inside mesoporous 

weakly depend on the vapour density. 

glasses 

The decrease of diameter with increasing velocity is 
physically expected, although measurements in the water- tsacKgrouna 

air system generally find a stronger V^"^ dependence '28]. rj^-, , rij.- j i.- c a 

■'° ■' b V f --"ihe phenomenon of condensation and evaporation of nu- 

However, this conclusion is at variance from the results •, . j- • -^i • r n- • ^i 

' ids m mesoporous media, i.e. with pores sizes lalhng m the 

of direct measurements of the diameter distribution on • ^ j- j 1 r 1 ^ • 

nm-/ini range, is studied as an example 01 a phase transi- 

the pipe axis using the PDPA [22ll23l . Thev show a nearlv ,■ ■ -cj- j j c j-mrn j j 

^ ^ ° ' — " — ' -^ ■' tion m presence 01 disorder and confanement [29j , and used 

exponential size distribution, with an average diameter , , , . j- j_ -u i- t-> ^1 4-1 j ^■ 

^ ' ° to characterise pore distributions. Both the condensation 

increasing with the vapour velocity (approximately from , , . ■ 1 ^i r ^- r j 

° ^ J V i-jr- J g^Q^j evaporation processes involve the formation of do- 

20 to 40 um). Measurements at temperatures between 1.8 ■ r i- • 1 • • i rm 

' ' ^ mams ot liquid or vapour on a microscopic scale, the as- 

and 2.0 K show that the diameter dependence on vapour ■ u- j a ^ j-- r ^-i-j i^ i-i^ i-i^ 

^ ^ sociated fiuctuations 01 optical index scatter light so that 

density follows a similar trend. Both points are inconsis- •,■,,,•,,1,11 • , 

■' ^ it IS attractive to study the phenomenon m transparent 

tent with the above analysis of the stratification. The rea- ,. , .,. , ,^r ■^■ 

porous media, such as silica glasses (Vycor, silica aero- 
son for this discrepancy remains to be understood. , n tt • • ^ ^i 1 ^i 1 r 
^ ■' gels, ... J. However, m some instances, the length scale of 

Despite this problem, the range of average diameters the fluctuations can reach several hundreds of nanome- 

given by the PDPA (« 10 — 30/im) is consistent with the ters, impeding single-scattering measurements for usual 

angular dependence of the intensity scattered at small an- fluids. This problem can be solved by using helium as a 

gle (figure[3]), and, independently, on the value of S, com- fluid. Figure [1] shows that, for helium droplets of diameter 

bined with the concentration of droplets, as inferred from 200 nm, the mean free path remains larger than one mil- 

the PDPA count rate and velocity. These diameters are limetre even for an extrapolated volume fraction of 50%, 
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allowing conventional light scattering measurements on a as long as its neighbours remain filled, so that desorption 

millimetre thick sample. This would not be possible for, should take place through a collective percolation process 

e.g., water, where the mean free path for a similar distri- starting from the surfaces of the sample. 

bution of droplets would be 200 times smaller. We have 

These different explanations for hysteresis predict dif- 

taken advantage of this specificity of helium to address 

ferent shapes for the hysteresis loop. For example, the per- 

several fundamental problems in the field using optical 

eolation mechanism should manifest through a sharp kink 

means. 

in the desorption isotherm at the percolation threshold|32|. 

The processes of condensation and evaporation are usu- followed by a nearly vertical portion, as observed in Vy- 
ally characterised by sorption isotherms, which measure cor and other materials. Also, for weak enough disorder, 
the adsorbed amount of fluid as a function of increasing t^e disorder-based mechanism predicts a change of shape, 
(for condensation) or decreasing (for evaporation) pres- ^0™ smooth to steep, of the adsorption isotherm as the 
sure. In the case where the dense phase of the fluid wets temperature is decreased below a critical value |33]. Still, 
the substrate, one observes the reversible adsorption of a measuring the shape of the loop only is ambiguous, and 
thin film at low vapour pressure, followed by capillary con- optical measurements allow to test more precisely the dif- 
densation, a rather abrupt filling at a pressure of order, ferent mechanismes proposed. This is beautifully iUus- 
but smaller, than the saturated vapour pressure. The later trated by the work of Page et al, who studied the conden- 
process is hysteretic, emptying occurring at a lower pres- sation and evaporation of hexanc in Vycor using light scat- 
sure than filling. A central question is the origin of this terin^, and found long-range correlations in the vapour 
hysteresis. Different explanations invoke the metastability distribution, characteristic of a percolation process, close 
of the gas-like phase[30|, changes in the shape of menisci to the kink of the desorption isotherm[34J. However, up to 
between condensation and evaporation 3 Ij, collective ef- now, there were no systematic study of the correlation be- 
fects for networks of connected pores such as Vycorp], tween the shape of the isotherms, and the distribution of 
and, more recently, the energetic and geometric disorder Auid inside a porous medium, both at the microscopic and 
of the porous media[35]. A characteristic feature of the macroscopic levels. This is the aim of our experiments. In 
two last explanations is that condensation and evapora- 



tion differ in a fundamental way, which is the infiuence of 2 ^^^^ ^.j^^^^ ^^ ^^^^ experiment, results were limited to pres- 
nucleation. Unlike condensation, which can proceed from gures very close to the percolation threshold. For smaller pres- 
the adsorbed film, evaporation requires the nucleation of sures, multiple scattering was so large that no light was trans- 
vapour in the dense phase, which involves an activation mitted, in agreement with the discussion at the beginning 
energy. If this process cannot occur, a pore cannot empty of ^3.11 
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Fig. 12. Illustration of the scattering geometry. The sample 
is a disk, or (here) a fragment of a disk, illuminated at 45° 
from its faces by a thin laser sheet. For these pictures, the 
background is illuminated with a weak extended laser beam 
through the viewport normally used for observation at 135°. 
(a) corresponds to the filled aerogel (the bulk interface is visible 
below the aerogel), and (b) to the beginning of the desorption 
process. In the filled state, scattering is only due to the silica 
backbone, immersed in liquid helium at an uniform density. 
The interception lines of the laser sheet with the back and 
front surfaces are brighter , due to forward scattering by these 
surfaces. During desorption, white regions appear, due to the 
coexistence of liquid and vapour on a microscopic scale. 

the following, we will illustrate the potential of our tech- 
nique by several examples. 

3.2 Experimental 

We study different porous glasses, silica aerogels and Vy- 
cor, using an optical cryostat with optical ports 45° apart. 
The sample cell, 20 mm in diameter and 4 mm thick, has 



two sapphire windows, allowing observation under differ- 
ent scattering angles. Sorption isotherms are measured be- 
tween 4.2 K and the critical bulk temperature (Ri5.2 K), as 
described in |36j and [37j . Simultaneously, the distribution 
of helium inside the sample is determined by light scatter- 
ing measurements. To that aim, the sample is illuminated 
by a thin He-Ne laser sheet under a 45° incidence with re- 
spect to its faces, and imaged at different angles (45°, 90°, 
and 135°) using CCD cameras, as illustrated on figure [T2l 

3.3 Evidence for a disorder-driven transition in silica 
aerogels 

Recent numerical studies |551l38) . based on a mean- field 
density functional theory, suggest that hysteresis could re- 
sult from the disorder of the porous media. They predict 
that a disorder-driven transition could occur as a function 
of disorder or temperature, similar to that occurring in 
the Random Field Ising Model at zero temperature |39j . 
This out-of-equilibrium phase transition implies a change 
of shape, from smooth to steep, of the condensation branch 
of the hysteresis loop, when the porosity is increased at a 
constant temperature or, alternatively, when the temper- 
ature is decreased below some critical value, at a constant 
porositv[iD]. 

We have recently provided the first evidence for this 
transition in a base-catalyzed silica aerogel[35i. In this 
sample (denominated BlOO), the silica forms a complex 
arrangement of interconnected strands, resulting in a large 
porosity (95%) and an associated weak disorder. For this 
porosity, the recent theoretical studies, performed on aero- 
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Fig. 13. Isotherms for a base-catalyzed aerogel of porosity 
95% (BlOO). The adsorbed fraction in aerogel, (p, is plotted 
as a function of the pressure, referred to the bulk saturated 
vapour pressure. The scanning time is typically one day for 
each branch of the loop. The adsorption branch evolves from 
steep to smooth as the temperature is increased, whereas, for 
desorption, it is nearly vertical at all temperatures. 




Fig. 14. Images of aerogel BlOO, observed at 45° from the 
incident laser sheet, for increasing adsorbed fractions #. The 
condensation morphology is correlated to the steepness of the 
adsorption isotherm (figure [13]). 

gels numerically synthesised by Diffusion Limited Clus- 
ter Aggregation [ID] . predict the occurrence of a disorder- 
driven transition. As shown by figure I13[ the adsorption 
isotherms become steeper at low temperatures, which is 



indeed in agreement with the scenario of a disorder-driven 
transition, and cannot be explained by the usual descrip- 
tion of the capillary condensation. The corresponding op- 
tical observations are shown in figure [T4| At low adsorbed 
fraction <P {<P < 50—60%), the scattered intensity increases 
uniformly, corresponding to the development of liquid do- 
mains on a microscopic scale (up to typically 200 nm, 
from the absolute value of the scattered intensity and its 
angular dependence). Above ^ = 50 — 60%, the behavior 
depends on temperature. For the low temperature points 
(corresponding to the steep isotherms), the injected he- 
lium completely fills some regions of the aerogel, making 
them dark. The size of these dark regions increases up to 
the point where the whole aerogel is filled. In contrast, 
at temperatures where the isotherm is smooth, the con- 
densation proceeds in an homogeneous way. As discussed 
in detail in ref. 35J , these observations are consistent with 
the occurrence of a disorder-driven phase transition. 

3.4 Absence of long-range correlations during 
desorption 

We have also studied the behavior during desorption[41j. 
In order to compare to the experiment of Page et al, 
we have studied a sample of Vycor. At the lowest tem- 
perature presently achievable with our set-up (4.3 K), 
the desorption isotherm does not have an initial verti- 
cal portion (figure llSp , in contrast with the situation at 
lower temperatures 42 . The scattered intensity decreases 
monotonously from the empty state to the filled state. 
This means that, unlike in aerogels, the initial adsorption 
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Fig. 15. Hysteresis loop for Vycor at 4.34 K. Note that hys- 
teresis occurs over a small range of <f . The inset shows the scat- 
tered intensity at 45° as a function of $. Along the hysteresis 
loop, scattering is slightly larger for adsorption (squares) than 



it is for desorption (circles). 
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Fig. 16. Images of desorption for aerogel BlOO, at 4.47 K and 
5.08 K, for decreasing liquid volume fractions $. Desorption 
starts from the aerogel surfaces. 

of the film decreases, instead of increases, the scattered in- 
tensity. This is related to the difference in porosities. For 
the aerogels, the silica concentration is small (several %). 
Adding more silica or a film of helium on the surface of 
the filaments increases the scattered field. In contrast, for 



Fig. 17. Hysteresis loop for aerogel B55 at 4.70 K and cor- 
responding images during the desorption process. As in BlOO, 
desorption starts abruptly and close to the surfaces. 

Vycor, the silica concentration is large (70%), so that one 
might see the scattering as originating from the cavities. 
Reducing the cavity size decreases the scattered signal, 
and so does the addition of an helium film. More sur- 
prizingiy, the scattered signal never exceeds its value in 
the empty state, even along the hysteresis loop where it 
is slightly larger for adsorption than for desorption (in- 
set of figure [T5)) . This implies that no microdomains are 
formed on the scale of the hundred of nanometers. More- 
over, there is no evidence for the long-range correlations 
during desorption which would sign a percolation process. 
On the other hand, the kink in the desorption isotherm 
is not a sharp one, so that our Vycor experiment does 
not rule out the usual hypothesis of a correlation between 
the desorption mechanism and the shape of the desorption 
isotherm. 

What about aerogels? In this case, the desorption iso- 
therms present a vertical portion for all temperatures (fig- 
ure [T5)) . The same holds for the two other aerogels de- 
scribed in ref.l35l. which have different microstructures 
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and porosities. Based on the above hypothesis, one would 
conclude that a percolation process is at play. Figure [TBI 
shows pictures of the desorption process for the previous 
sample. When the adsorbed fraction starts to decrease 
(corresponding to the kink in the isotherms), bright re- 
gions appear, demonstrating the apparition of vapour "bub- 
bles" at the microscopic scale. In the initial stage of emp- 
tying, these regions are located close to the surfaces of 
the aerogel, (i.e. near the bright lines marking the back 
and front surfaces, and on the bottom part of the aero- 
gel). The same conclusions apply to the other aerogels. 
For example, figure [17] shows the desorption process for 
a base-catalyzed sample of porosity 97.5% (B55). Desorp- 
tion initially occurs along some "fingers" , some of which 
distinctly originate from the aerogel surfaces. This points 
to a specific role of surface during the desorption process, 
seemingly consistent with the percolation scenario. 

In this scenario, one expects a large increase of scatter- 
ing with respect to the case of condensation, as observed 
for hexanc in Vycor[34j. However, for all three samples 
studied, we found no evidence for such an increase. This 
is illustrated by figures [H] and ^^ Figure [TH] shows, for 
the base-catalyzed sample BlOO, the intensity scattered 
at 45° from a point close to the surface, where the bright 
regions first appear during desorption. At a given average 
adsorbed fraction, the signal during evaporation is some- 
what larger than during condensation, but by a factor less 
than 3, which can be explained [41] by a modest increase in 
size (typically 40%) of the correlated domains. The case of 
the neutrally-catalysed sample (N102) is even more spec- 
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Fig. 18. Intensity scattered at 45° for BlOO. The intensity 
is the average grey level in the rectangle of the picture, nor- 
malized by the silica background in the same rectangle. The 
difference between adsorption and desorption is moderate, im- 
plying that the correlations during desorption do not extend 
on much larger scales than during adsorption. 



tacular. Despite a sharp kink in the desorption isotherm, 
the optical signal shown in figure [19] is nearly reversible as 
a function of <?. The same results hold for all temperatures, 
and for B55, which means that, for all three aerogels, des- 
orption does not involve much longer length scales than 
adsorptior[f]. To our knowledge, these results are the first 
evidence that the kink in the desorption isotherm docs not 
necessarily follow from a percolation process. 



■^ This contrasts with the results of neutron scattering studies 
for aqueous gelsi43j. However, the length scales probed in that 
case did not exceed a few tens of nm, an order of magnitude 
smaller than those we probe. 



P.E. Wolf et at Probing helium interfaces with light scattering 



19 




Fig. 20. Evidence for bulk nucleation during desorption in aerogel B55: (a) Several slices at constant liquid volume fraction 
and temperature; (b) Reconstructed volume seen parallel to the surface of the disk. The circled regions do not contact the outer 
perimeter of the disk; (c) View orthogonal to (b), showing that the circled regions do not contact the disk faces either (top and 
bottom on (c)). (b) and (c) are reconstructed using the ImageJ plug-in Volume Viewer. 



3.5 Testing the occurrence of nucleation 

The absence of long-range correlations for BlOO contrasts 
with the predictions of the mean-field density functional 
approach. Numerical calculations of the structure factor 
for a base-catalysed aerogel of porosity 87% predict that 
the scattering during desorption should be much larger 
than during adsorption [33]. A possible explanation would 
be that, in the real aerogel, thermal activation, which is 
neglected in the numerical studies, could allow to over- 
come the nucleation barrier. If this were true, one should 
observe the apparition of vapour "bubbles" , hence of bright 
regions, in the bulk of the sample. For testing this issue, 
a single slice of the sample is not enough. For example, 
the bright "fingers" which appear in the middle of the 
laser sheet in figure [17] could be connected to the surfaces 
outside of the sheet plane. Detecting the occurrence of nu- 
cleation thus requires to reconstruct the whole 3D sample. 
This is possible by scanning the laser sheet parallel to it- 



self and capturing the resulting images. In this way, one 
obtains a stack of parallel slices of the sample, which can 
be reconstructed using the freeware ImageJ |3S]. Figure!^ 
shows B55, soon after the beginning of desorption. Several 
small bright regions can be identified as non-connected 
to the surfaces. If confirmed, this result will be an evi- 
dence for nucleation. Light scattering could then be used 
to study how the nucleation mechanism depends on the 
temperature and the microstructure of samples. 

4 Conclusions 

In this paper, we have shown how the weak optical con- 
trast of helium droplets or bubbles allows to perform quan- 
titative light scattering measurements on diphasic systems 
at large volumic fractions. We have exploited this property 
to study optically the atomisation of liquid helium and the 
gas-liquid phase transition in porous media. The original 
results obtained so far illustrate the richness of blending 
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Fig. 19. Intensity scattered at 135° for N102. Due to its larger 
correlation length, the empty aerogel N102 scatters light much 
more strongly than BlOO or B55. The scattering mean free 
path is comparable to the sample thickness, so that the light 
backscattered at 135° is significantly attenuated, depending on 
the probed depth. To avoid the complication of this effect, the 
scattered intensity is measured close to the entrance surface of 
the beam. When plotted as a function of $, it is nearly iden- 
tical between adsorption and desorption. This is valid as well 
for other points in the sample, showing that, on the scale of 
the hundred of nanometers, the distribution of liquid is unex- 
pectedly the same for adsorption and desorption. 

low temperature and soft matter physics, and will be a 
motivation to pursue this approach further. 
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